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Abstract 



We present predictions for the diffractive production of heavy quarks in deep- 
inelastic scattering and hadron-hadron collisions and for double diffractive 
dijet production in hadron-hadron collisions. With the assumption of hard 
scattering factorization, the predictions are made using diffractive parton den- 
sities that we have previously fitted to HERA data. Comparisons of our pre- 
dictions with preliminary data from HERA and the Tevatron indicate that 
factorization is obeyed in diffractive DIS but fails in hadron-induced processes. 
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I. INTRODUCTION 



In a previous paper [1], we obtained diffractive parton densities by fitting diffractive deep 
inelastic scattering (DIS) and photoproduction data from HERA, and then we used these to 
predict diffractive production of 1^'s and dijets at the Tevatron. The underlying principle 
behind these calcTiIations was the notion of hard scattering factorization, a property which 
we aimed to test. 

A factorization theorem has been proven [2] for diffractive DIS and direct photoproduc- 
tion processes. The proof establishes the universahty of diffractive parton densities for the 
class of processes to which the theorem applies. This means that one may extract the parton 
densities from a subset of data, and then, using these same densities, reliably predict other 
diffractive DIS and direct photoproduction distributions. The results of our work in [1] are 
consistent with this theorem. 

In contrast, no factorization theorem has been established for diffractive hadron-hadron 
processes, and in fact, there have been several arguments put forward in favor of nonfac- 
torization [3,4] in this case. Again, our results in [1] were consistent with this, when we 
compared our predictions for hadron-hadron processes to data. 

In this paper, we seek to further extend our tests of factorization. Using the same diffrac- 
tive parton densities that we extracted in [1], we present calculations of cross sections for 
one process where factorization is predicted to be valid — diffractive heavy quark production 
from DIS — and for two processes where factorization is predicted to fail — diffractive heavy 
quark production in pp interactions, and diffractive dijet production via double Pomeron 
exchange. We will compare our predictions with preliminary data from HERA [5,6] and the 
Tevatron [7,8]. 

Apart from the general desire to test hard-scattering factorization, a particular motiva- 
tion for studying diffractive DIS charm production is that we expect charm production to 
be a substantial fraction of the diffractive DIS cross section. (Recall that in inclusive DIS, 
the charm structure function F^"^ is ~ 25% of the inclusive F2.) Moreover, since open charm 
production is dominated by the photon-gluon fusion process, it also provides a good test of 
the normalization of the diffractive gluon density. A particular result of our work in [1] was 
that the diffractive ghion density is much larger than the diffractive quark density. This was 
most strongly driven by the photoproduction data. The hadron-induced processes that we 
study in this paper are also highly sensitive to the diffractive gluon density. 

The paper is organized as follows. In section II, we discuss the diffractive parton densities 
and cross section formulae that we will use in the calculations. Our numerical results are then 
presented and discussed in relation to preliminary data in section HI. We then summarize 
our findings and give our conclusions in section IV. 

II. CROSS SECTIONS 

Since our goal is to test hard scattering factorization, we will use, in our calculations 
of the diffractive cross sections, the parton densities which we previously obtained [1] from 
fits to diffractive DIS and photoproduction data from ZEUS and HI. Similar fits were also 
obtained independently in [9,10], using only data from ZEUS. Note that there have also 
been previous fits to diffractive DIS data alone [11]. 
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In [1], where we assumed that the diffractive parton densities can be adequately repre- 
sented as a Pomeron flux factor times parton densities in the Pomeron, the flts are labeled 
A, B, C, D and SG. In fits A, B, C and D, the parton densities in the Pomeron all have the 
form, 

Pf,/AP, Qo) = a, - /3) + a, (1 - Pf] , 

/9/,/F(/9,Qo) = a,/3(l-/?), (1) 

at an initial scale Qo — 2 GeV, but with different constraints on the parameters. For 
example, in parameterizations A and C we required that there be no gluons in the initial 
parton densities, i.e., = 0. The super-high-glue fit, SG, has the same form of quark 
density as in Eq. (1) but with a gluon density that is peaked near /3 = 1 and is given by 

(3fM/3,Qo)-a,(3'(l-(3f\ (2) 

The appropriate Pomeron flux factor to use with the above parton densities is one due 
to Donnachie and Landshoff [12]: 

/p°/p(^p) = /° dt 

where rrip is the proton mass, Po — 1-8 GeV""*^ is the Pomeron-quark coupling and a{t) — 
ap + 0.25t is the Pomeron trajectory. 

The explicit parameters for our flts are presented in Table I. Note that flts A and C 
have low gluon content (i.e., gluons are only produced via evolution) while flts B, D and 
SG contain large amounts of initial gluons. In [1], we found that only the high-glue flts are 
able to flt the photoproduction data while the low-glue flts badly underestimate the same 
data. All flve parameterizations are in reasonable agreement with the DIS data, at least as 
regards normalization. 

Our motivation for working with flts that do not in fact flt all the data is that the gluon 
distributions that are needed to provide good flts to photoproduction data are an order of 
magnitude above the quark distributions. This is quite unlike the usual inclusive parton 
densities in the proton. Thus, a comparison of results obtained with the low-glue and high- 
glue flts should dramatically distinguish those processes that are particularly sensitive to 
the diffractive gluon density from those that are merely sensitive to the gluon densities in 
their flner details. 



4m^ 



2.8t 



Aml-t \l-t/Q.7 



X 



l-2a(t) 



(3) 



A. Heavy quarks in DIS 

We compute the total cross section for diffractive production of heavy quarks in DIS 
with the usual formula 

= / dXadXt>fa/F{Xa,n)fF/pixF)^ar^ (4) 

where Xa and xp are momentum fractions of parton a and the Pomeron, respectively, 
/a/p(a;o,/x) is the distribution function at scale of parton a in the Pomeron, f-p/p{x-p) 
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Fit 


ttq 




dq 


A 


0.240 ± 0.006 








B 


0.239 ±0.006 


4.5 ± 0.5 





C 


0.249 ±0.011 





-0.031 ± 0.029 


D 


0.292 ±0.013 


9.7 ± 1.7 


-0.159 ±0.029 


SG 


0.225 ± 0.008 


7.4 ± 2.2 






TABLE I. Fit parameters and associated errors for the diffractive parton densities with ap = 1.14 
as presented in [1]. 



is the Pomeron flux factor and da^y* is the coefficient function for parton- virtual photon (7*) 
scattering. The next-to-leading order form of this function may be found in [13]. The sum 
in Eq. (4) runs over all active parton flavors. 

The calculation that wc have used is based on HVQDISvl.l by B. W. Harris and J. Smith 
[14] which computes inclusive (i.e., without the diffractive requirement), fully differential 
heavy quark distributions in DIS both to leading (LO) and next-to-leading order (NLO). A 
Peterson fragmentation routine is built-in and can be applied to the flnal state heavy quark 
four-vectors to simulate hadronization effects. We modifled this calculation to compute 
diffractive cross sections using the diffractive parton densities obtained in [1] . 

B. Heavy quarks in hadron collisions 

The cross sections for heavy quark pair (QQ) production in pp collisions for the case 
when the p diffracts are calculated using 

^QQ = ^ J dXb dXp h/p{xp) fb/FiXb, fa/pi^a, ^ab^Qqi^aXbS, UIq, /X^), (5) 
a,b 

where fa/p{xa,fJ^) is the proton parton distribution function, s is the center of mass energy 
squared, uiq is the heavy quark mass and Ofja^QQ is the partonic hard scattering function 
for heavy quark pair production. The analytic expressions for this function up to NLO are 
given in [15]. 

It is useful to present the above cross section as a fraction of the inclusive cross section. 
Thus, we will also compute the diffractive rate R — crQQ/<^QQ-, where 
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C. Dijets via double Pomeron exchange 



For the production of dijets via double Pomeron exchange in pp colhsions, we compute 
the cross section by using 

^DPE ^ ^ j dXa dXb dXw/p dXw/p fw/p{xp/p) fa/w{Xa, jj) f¥/p{xw/p) fb/w{Xb, jj) ^ab, (7) 
a,b 

where Uab is the partonic 2 — > 2 scattering function. The LO expressions for Uab may be 
found in [16]. 

We will also compute the ratio of the above double diffractive cross section to that for 
single diffraction {(t^^). The latter, assuming the p diffracts, is calculated with the formula 




III. RESULTS AND COMPARISONS WITH DATA 

In computing the diffractive cross sections, we used the diffractive parton densities dis- 
cussed in section II. These are evolved using full NLO QCD, via the CTEQ package [17], 
from the initial scale Qq to the appropriate scale in the particular calculation. 

For the calculations involving heavy quarks, NLO hard scattering matrix elements [15] 
were used together with a two-loop running coupling a^. To be consistent with the rcnor- 
malization scheme used in [15], one needs to use the coupling which is valid for three and 
four flavors for charm (c) and bottom (6) production, respectively. At NLO, the matching 
condition between the 3- and 4-flavor couplings at = rric is that the coupling is continu- 
ous, and similarly for the transition from 4 to 5 active flavors — see [18]. For consistency in 
the calculations of 0"^^', we also used proton parton densities with the appropriate number 
of flavors [19] in the flxed flavor scheme [20]: CTEQ4F3 (CTEQ4F4) for charm (bottom) 
production. 

A. Diffractive charm production at HERA 

Both the HI [5] and ZEUS [6] collaborations have reported preliminary data on diffractive 
charm production in deep inelastic scattering of 27.5 GeV positrons and 820 GeV protons. 
They studied charm quark production, by tagging D* mesons through their decay, via the 
process 

e+ + p ^ e+ + p + {D*^ (D° X-7r+)7r±) + X. (9) 
1. HI data 

HI has measured a preliminary cross section for diffractive D* production shown in the 
last column of Table II. This is subject to the following cuts: 10 GeV^ < < 100 GeV^, 
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A 


B 


C 


D 


SG 


HI data 


a 


4.22 


196 


4.15 


419 


70.3 


380 It^lj (stat.) (syst.) 



TABLE II. Diffractive D* cross sections (pb) in DIS using HI cuts. The last column shows the 
preliminary measurement by HI [5]. 





xf^ = 0.01 


x^^ = 0.05 


xf^ = 0.1 


A 


0.74 


5.15 


8.67 


B 


76.7 


242 


315 


C 


0.75 


5.09 


8.46 


D 


165 


518 


670 


SG 


43.2 


87.4 


105 


ZEUS data 


604 ±171 (stat.) Ilifs (syst.) 



TABLE III. Diffractive D* cross sections (pb) in DIS using ZEUS cuts and for several different 
limits on Xf>. The bottom row shows the preliminary measurement by ZEUS [6]. 

Pt{D*^) > 1 GeV, 0.06 <y < 0.6, xr < 0.05 and \viD*^)\ < 1-25, where is the photon 
virtuahty, y is the fractional energy loss of the incident e"*" in the proton rest frame while 
Pt and t] are the transverse momentum and pseudorapidity, respectively, of the D* meson. 
These data include proton dissociation with masses up to 1.6 GcV. 

The NLO cross sections that we calculated with the same cuts are also given in Table 
II. These values were obtained with Peterson fragmentation switched on in the code with 
Ec — 0.035 and the probability for the charm to fragment into a D* meson, P(c —>■ D*), 
set equal to 0.26. The renormalization and factorization scales in Eq. (4) were both set to 
/X = ^/Q"^ + 4m^, where the mass of the charm quark, nic, is set equal to 1.5 GeV. 

We see that with the low-glue fits A and C, our predicted cross sections are two orders of 
magnitude smaller than the data. With fits B and SG, our cross sections are about two and 
five times smaller than the central data value but within the errors. Fit D, which provided 
our best fit to the DIS and photoproduction cross sections, yields a diffractive charm cross 
section which is in best agreement with the data, being within 10% of the central value. 

2. ZEUS data 

For the same process, the cross section obtained by ZEUS, corrected for a (31 ± 13)% 
double dissociation contribution, is given in the bottom row of Table III. The following 
experimental cuts apply for this measurement: 10 GeV^ < < 80 GeV^, Pt{D*^) > 
1 GeV, 0.04 <y< 0.7, \viD*^)\ < 1.5, and ^^.^ax < 2. 

Using the same cuts, we derive the cross sections shown in Table 111. We used the same 
fragmentation parameters and choice of scale as in the calculations for comparison with the 
HI data. Since the ZEUS data have an rj^ax cut, it is not clear what limit to impose on xp. 
Hence, we have chosen several different values, as shown. 
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We find that if wc use tlie cut < 0.01, tlie predictions badly underestimate tlie data. 
Using cuts of 0.05 and 0.1, the predictions obtained with the low-glue fits A and C are still 
one to two orders of magnitude smaller than the data. The predictions using the high-glue 
fits B, D and SG have better agreement, with those obtained from D being closest to the 
data. 

3. Sizes of predictions 

We make several remarks on the relative sizes of our predictions for this process. The 
cross sections derived with the low-glue fits (A and C) arc from 16 to about 100 times smaller 
than those obtained from the high-glue fits B, D and SG. This large discrepancy is due to 
the fact that the production of heavy quarks in DIS is dominated by the "y* — g channel, 
with quark contributions entering only at next-to-lcading order. With fits A and C, gluons 
are generated only via evolution of the parton densities. Therefore, the predictions using 
fits A or C arc suppressed by at least a factor of relative to those with the high-glue fits. 
Furthermore, another factor of 10 suppression results from the quark coefficients in Table I 
being an order of magnitude smaller than the gluon coefficients. 

Among the high-glue fits, the cross section using SG is smallest. This indicates that, 
for this process, the bulk of the contribution comes from partons with momentum fractions 
(x) that are moderately far from unity; the gluon density of SG dominates that of B or D 
only in the region a; ~ 1. The cross sections using fits B and D have about a factor of two 
difference. This is due to the gluon coefficient ag for D being about twice that for B and the 
dominance of gluon-induced scattering in this process. 

4- Differential cross sections 

Figure 1 shows our predicted differential cross sections using HI cuts plotted as functions 
of W (the center-of-mass energy of the virtual photon-proton system), Q"^, pt{D*), r]{D*) 
and xp. The same comments apply regarding the relative sizes of the predictions for each 
distribution as those made above. The plots strongly illustrate how well diffractive charm 
production in DIS can determine the gluon density as evidenced by the large difference 
between the predictions using the realistic high-glue fits (D, B and SG) and the low-glue 
ones (A and C). 

B. Heavy quark production from hadrons 

At this time, there are no data available for comparison with our predictions although 
there is a preliminary report from CDF [21] which quotes that "(0.18 ± 0.03)% of central 
high— b/c are diffractively produced" in pp interactions at ^/s = 1800 GcV. One may 
thus expect that data on diffractive heavy quark production will soon become available. It 
is useful, therefore, to make some predictions to study the differences between the diffractive 
parton densities. 

For open heavy fiavor (QQ) production from pp interactions at y/s = 1.8 TeV, we 
obtained the NLO diffractive cross sections (assuming that only the p diffracts) shown in 
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50 100 150 200 250 ^ 

W (GeV) Q' (GeV') 




FIG. 1. Various distributions for diffractive D* production in DIS using HI cuts. Solid curves are 
for predictions obtained with fit D, dotted for those with fit B, heavy dashed for fit SG, dot-dashed for 

fit C and dashed for fit A. (The last two lie almost on top of each other.) Only fits B,D and SG are to 
be regarded as realistic predictions. The two plots at the bottom are identical except for the choice of 
scale for the axes. 
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^max 
X jp 


= 0.01 




= 0.05 


^max 


= 0.1 






R 




R 




R 




A 


5.2 


0.2% 


8.6 


0.4% 


10.5 


0.5% 




B 


194 


8.8% 


323 


14.6% 


390 


17.7% 


Q — charm 


C 


5.0 


0.2% 


8.1 


0.4% 


9.9 


0.5% 




D 


414 


18.7% 


689 


31.2% 


831 


37.6% 




SG 


69.1 


3.1% 


108 


4.9% 


128 


5.8% 




A 


0.1 


0.1% 


0.2 


0.3% 


0.3 


0.5% 




B 


1.7 


2.6% 


4.5 


6.9% 


6.4 


9.9% 


Q — bottom 


C 


0.1 


0.1% 


0.2 


0.3% 


0.3 


0.5% 




D 


3.5 


5.4% 


9.5 


14.6% 


13.5 


20.8% 




SG 


0.8 


1.2% 


1.9 


2.9% 


2.6 


3.9% 



TABLE IV. Diffractive cross sections crgg (/xb) and rates R = (^qq/(^qq for open charm and 
bottom production at the Tevatron using several different cuts on xp. These values are valid for the 
case when only the p diffracts. 

Table IV. The corresponding rates, R = <^qq/<^qq, are also shown in the same table. We 
used rric = 1.5 GeV and rrib = 4.5 GeV for the charm and bottom quark masses, respectively, 
and set the scale /j, in Eqs. (5,6) to the heavy quark mass. 

As in the case of heavy quark production in DIS, the predicted cross sections using the 
(unrealistic) low-glue fits A and G are one to two orders of magnitude smaller, depending on 
the cut on xp used, than those obtained with fits B, D and SG. Again, this is because most 
of the cross section for this process results from gluon-induced {g — g and g — q) scattering. 

Prom Table IV, we see that for charm (bottom) production, the predicted diffractive 
cross sections using the high-glue fits are about 3% to 40% (1% to 20%) of the inclusive 
cross section. The low-glue fits result in much smaller rates: 0.2% to 0.5% and 0.1% to 0.5% 
for charm and bottom production, respectively. 

C. Dijets via double Pomeron exchange at the Tevatron 

CDF has preliminary data [7,8] on double diffractive dijet production frompp interactions 
at y/s — 1.8 TeV which are characterized by central dijets bounded by rapidity gaps near 
the direction of each incoming hadron. They used a Roman pot to trigger on the diffracting 
antiproton and tagged on rapidity gaps along the direction of the diffracting proton. The 
measured cross section [7] for dijets with transverse energy Ex > 7 GeV and cuts of 0.05 < 
xp/p < 0.1 and 0.015 < x^/p < 0.035, where x^/h is the fractional momentum lost by the 
diffracting hadron h, is shown in the last column of Table V. 

There also exist data on the double diffractive cross section cr^^^ expressed as a fraction 
of the single diffractive cross section a^^. The fraction cr^^^/a^^ measured by GDF [8] is 
shown in Table V. In this case, the denominator is the cross section for production of 
dijets with Ex > 7 GeV and a cut 0.05 < x^/p < 0.1 (obtained with a Roman pot trigger 
on the diffracting antiproton). 
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A 


B 


C 


D 


SG 


Preliminary CDF data 




10.1 


898 


9.82 


3713 


206 


13.6 ±2.8 (stat.) ±2.0 (syst.) 




0.18% 


1.66% 


0.19% 


3.42% 


0.89% 


[0.17 ± 0.036 (stat.) ± 0.024 (syst.)]% 



TABLE V. Dijet cross sections in pp interactions via double Pomeron exchange using CDF cuts. 
The first row shows the cross section (nb) while the second row presents the rates expressed as a fraction 
of the single diffractive cross section (o"^^). The last column shows the CDF preliminary measurements. 



The cross sections cr^^^ and ratios a^^^/a^^ that we obtain using the same cuts and 
the diffractive parton densities shown in section II are also presented in Table V. In our 
calculations, wc set the scale fi in Eqs. (7,8) equal to Et- 

Compared to the data, we find that the predicted cross section a^^^ and ratio cr^^^/ a^^ 
obtained with the low-glue fits, A and C, agree well with the data. However, as we have 
seen, these fits substantially underestimate diffractive photoproduction and diffractive charm 
production. With the reahstic high-glue fits B, D and SG, the predictions badly overestimate 
the data, being one to two orders of magnitude larger. 

The relative sizes of the predicted cross sections in Table V can be explained in the 
following way. As in the case of heavy quark production, dijet production is dominated by 
gluon-induced scattering at the parton level. However, in this case, quarks do contribute at 
the leading order. Even then, the cross sections using fits A and C are much smaller relative 
to those obtained with the fits B, D and SG than in the heavy quark case (see Tables II, 
III and IV). The reason is that in the double Pomeron exchange cross section, there is an 
extra factor of the gluon density for the second Pomeron. Fits A and C are then expected to 
produce smaller cross sections because they can only supply gluons for the hard scattering 
via evolution. Similarly, the factor of four difference between the results for fits B and D is 
a consequence of fit D having about twice as much initial gluon as B. 

IV. CONCLUSIONS 

We have presented diffractive heavy quark and dijet cross sections calculated with the 
assumption of hard scattering factorization. Diffractive parton densities that were fitted [1] 
to DIS and photoproduction data from HERA were used in these calculations. In [1] , our 
best fit to the data resulted in fit D. 

We find that in the case of diffractive D* production, the predicted cross sections using 
the best-fit parton densities, fit D, also have the best agreement with the data. This strongly 
supports the notion of hard scattering factorization for diffractive lepton-hadron processes. 

For heavy quark production from pp collisions and assuming only the p diffracts, we 
obtained rates R = (^qq/cTqq, using the high-glue fits, from 3% to 40% and 1% to 20% for 
open charm and bottom production, respectively. The rates obtained with the low-glue fits 
are less than 1% for either charm or bottom production. It will be interesting to see how 
well these figures compare when data from the Tevatron become available. 

Our results for diffractive dijet production at the Tevatron indicate a breakdown of fac- 
torization for diffractive hadron-hadron processes. The predicted double Pomeron exchange 
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cross sections a with the high-glue fits are factors of 10-100 times larger than data. Sim- 
ilar results are obtained when comparing the ratio cr^^^/a^^ with data. Only the low-glue 
fits, which do not fit diffractive photoproduction data at all, yield numbers that agree well 
with the data. 
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